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ABSTRACT 

We use models of coupled thermal evolution and photo-evaporative mass loss to understand the formation 
and evo lution of the Kepler-36 system. We show that the large contrast in mean planetary density observed by 
ICarter e t al. (2012) can be explained as a natural consequence of photo-evaporation from planets that formed 
with similar initial compositions. However, rather than being due to differences in XUV irradiation between the 
planets, we find that this contrast is due to the difference in the masses of the planets' rock/iron cores and the 
impact that this has on mass loss evolution. We explore in detail how our coupled models depend on irradiation, 
mass, age, composition, and the efficiency of mass loss. Based on fits to large numbers of coupled evolution 
and mass loss runs, we provide analytic fits to understand threshold XUV fluxes for significant atmospheric 
loss, as a function of core mass and mass loss efficiency. Finally we discuss these results in the context of 
recent studies of the radius distribution of Kepler candidates. Using our parameter study, we make testable 
predictions for the frequency of sub-Neptune sized planets. We show that L8-4.0 planets should become 
significantly less common on orbits within 10 days and discuss the possibility of a narrow "occurrence valley" 
in the radius-flux distribution. Moreover, we describe how photo-evaporation provides a natural explanation for 
the recent observations of .Ciardi et al..(,2013J that inner planets are preferentially smaller within the systems. 



1. INTRODUCTION 



The Kepler-36 system (ICarter et al.ll2012h is a fascinating 
case-study for theories of planet formation and evolution. It 
contains a closely spaced pair of super-Earth mass planets 
with periods of 13.8 and 16.2 days orbiting a slightly evolved 
6.8 ±1.0 Gyr old Gl sub-giant that is 2.9 times more lumi- 
nous than the Sun. Although too faint for reliable radial veloc- 
ity measurements, the system exhibits strong Transit Timing 
Variations (TTVs), which allowed the pl anet densities to b e 
determined to better than 10% precision (ICarter et al.ll2012l) . 
Surprisingly, despite their extremely similar orbits, the plan- 
ets have densities that differ by almost an order of magnitude. 
The inner planet Kepler-36b has a mass of 4.5 ±0.3 and 
a density of 7.2 ±0.7 gcm~^, fully consistent with an Earth- 
like composition. Meanwhile the outer planet Kepler-36c 
has a mass of 8.1±g:^ but a density of only 0.87 ±0.05 
cm~^, meaning th at it must have a substantial H/He envelope 
ICarter etal.ll2012h . 

This poses an interesting problem for theories of planet for- 
mation and evolution: how is it that two planets in the same 
system with very similar highly irradiated orbits arrived at 
such radically different densities today? One possibility is 
that the planets simply formed in very different environments. 
Models of core accretion show that it is much easier to accrete 
a substa ntial H/He envelope w hen the local disk temperature 
is lower dlkoma & Horill20T2l) . Perhaps Kepler-36b formed at 
or close to its current orbit while Kepler-36c formed substan- 
tially further out and migrated inwards (Ida & Lin 2010). 

However, another possibility is that the planets did form 
in similar environments, but that subsequent evolution has 
caused them to diverge. In particular, photo-evaporation 
due to extreme ultra-violet (XUV) heating can remove large 
amounts of hydrogen/helium from highly irradiated planets 
through hydrodynamic mass loss. Models of XUV-driven 
mass l oss were first develop ed to study water loss from early 
Venus (iHuntenll982i:lKasting & Pollaca 19831) . and hydrogen 



loss f rom the early Earth (ISekiva et al.l Il980l IWatson et al.l 
119811) . These kinds of models wer e further developed to study 
mass loss from hot Jupiters (e .g., iLammer et al.l 120031 lYellel 
l2004tlMurray-Clay et al.ll2009l) . where there is clear evidence 
that atmospheric escape is an important physical process. 
This includes both trends in the population of h o t Jupiters 
dLecave her Des Etangsl 120071: iHu bbard et al.l l2 ()07blal: 
|Davis_& Wheatlejj |2009|; iHirenreich & Desert 2011] ) and 
direct detections of atmospheric escape ( Vidal-Madiar et aLl 
2003, 2004; Linskv et al. 2010; LecaveUer Des Etangs et M 
2010; Lecavelier des Etangs e t al.1120121: IHaswell et all 12011 
Ehrenreich et al. 2012). More recently, mass loss models 



have been used to study the new populations of super-Earths 
and sub-Neptunes being found by transiting surveys like 
CoRoT and Keple r (Jackson et al. 2012; Owen & Jacksoi 



■Ml i Lopez et al.ll2012t iWu & Lithwickll2012t lOwen & Wi 
2013). 
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Recently, in ILopez et alJ (120121) we showed that exoplanet 
compositions are subject to a photo-evaporation threshold. 
Observationally there are no planets with low bulk density 
and high incident flux, implying that low-mass planets with 
substantial H/He envelope s do not exist in this area of pa- 
rameter space. Following Lecavelier Des Etangs (2007) who 
studied hot Jupiters, Lopez et al. (2012.) were able to show 
that this threshold can naturally be ex plained as a critical 
mass loss timescale. Detailed models in ILopez et al.l (120121) . 
which include coupled thermal and mass loss evolution, re- 
produced this threshold with standard mass loss efficiencies. 
We further showed that many of the super-Earths and Nep- 
tu nes found by the Kep ler mission including those in Kepler- 
1 I dLissauer et al.ir201 ll) and Kepler-36 lie along this thresh- 
old, indicating that these planets may have undergone sub- 
stantial mass loss in the past. Here we show that photo- 
evaporation can be controlled by the mass of a planet's 
rock/iron core and that this provides a natural explanation for 
the divergent densities of Kepler-36 b&c. This process allows 
both planets to form with similar compositions and in similar 
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environments before being sculpted by their different mass 
loss histories. 

In addition to detailed studies of individual systems like 
Kepler-36, there is a growing body of literature examining 
the overall distribution of Kepler planet candidates. Detailed 
stu dies of planet occur rence rates by Petigura et al. (2013) 
and iFressin et al.l (120 13h have recently shown that is a sharp 
drop off in the frequency of planets above >2.8 7?®, at least 
within 50 and 85 days respectively. In contrast, planets 
with radii <2.8 Rof, seem to be equally common. Likewise, 
IHoward et al. (2012) and Youdin (2011) showed that there is 
a drop off in the frequency of candidates at extremely short 
periods <10 days; moreover, this drop off is particularly 
acute for su b-Neptune sized planets. In multi-planet systems, 
ICiardi et"aLl (i2013i) showed that there is tendency for outer 
planets to be larger than inner planets in the same system, 
particularly when those planets are smaller than Neptune and 
on orbits <20 days. 

Models have shown that photo-evaporation and thermal 
evolution can significantly alter the H/He inventories of super- 
Earth and s ub-Neptune sized planets (iBaraffe et al.l 120061; 
iLopez et al.i l2012t lOwen & Wul 120 13h . Moreover, planet 
structure models have shown that changing the mass of a 
planet 's H/He envelope has a large impact on it s resulting ra- 
dius ( Rogers & Seageril2010t ILopez et al.ll2012l) . As a result, 
models of thermal and mass loss evolution can shed light on 
many of the observed trends in the radius distribution of Ke- 
pler planet candidates. 

2. OUR MODEL 

For this work, we have used the co upled therma l evolu - 
tion and mass loss model presented in ILopez et"an (120121) . 
where additional model details can be found. Similar mod- 
els have been used to track the coupled evolution of rocky 
super-Earths (e.g , Jjac kson et al. 2010; Valencia et al. 2010|; 
Nettelmann et alJ 1201 ll) . hot N eptunes (e.g , Baraffe et al. 
20061). and hot Jup iters (e.g, IBaraffe et al., ,2004. ,2005; 



Hubbard etal.]l2007bllah . Be ginning shortly after the end of 



planet formation, we track planetary mass and radius as a 
function of age. The use of coupled model is essential, be- 
cause planetary radii are largest at young ages, when stellar 
XUV fluxes are highest. 

At a given age, a model is defined by the mass of its 
heavy element core, the mass of its H/He envelope, the 
amount of incident radiation it receives, and the internal en- 
tropy of its H/He envelope. Here we assume an isothermal 
rock/iron core with an Earth-like 2:1 rock/iron ratio, using 
the ANEOS olivine (Thompson 1990) and SESAME 2140 
Fe jLvon & Johnsonl ll992) equations of state (EOS). For the 
H/He envelope we a ssume a fully adiabatic interior using the 
iSaumon et alJ (119951) EOS. 

In order to quantitatively evaluate the cooling and contrac- 
tion of the H/He envelope, we use a model atmosphere grid 
over a range of surface gravities and intrinsic fluxes. These 
grids relate the surface gravity and internal specific entropy 
to the intrinsic flux emitted for a given model. These radia- 
tive transfer models are computed at a Uranus and Neptune- 
like 50 X solar metallicity atmosphere using the methods de- 
scribed in Fortney et al. (2007) and Nettelmann et al. (201 1 ). 
These atmosphere models are fully non-gray, i.e. wavelength 
dependent radiative transfer is performed rather than simply 
assuming a single infrared opacity. In addition, we include 
heating from radioactive decay in the rock/iron core and the 
delay in cooling due to the core's heat capacity. In order 



to correctly determine a planet's mass loss history, it is vi- 
tal to include these thermal evolution effects, since these will 
strongly affect a planet's radius over time. Radius, in turn, has 
a large impact on the mass loss rate as seen in Equation[T] 

Close-in planets like those in Kepler-36 are highly irradi- 
ated by extreme ultraviolet (EUV) and X-ray photons. These 
photons photo-ionize atomic hydrogen high in a pla net's at- 
mosp here, which in turn produces significant heating (i HuntenI 
119821) . If this heating is large enough, it can generate a hydro- 
dynamic wind that is capable of removing significant mass. 
We couple this XUV-driven mass loss to ourtherrnal e volution 
mode ls using the energy-limited approximation (iV^tson et alJ 
11981b . This allows a relatively simple analytic description of 
mass loss rates. 
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^t.de = (l-^ + ^) 



(1) 



(2) 



(3) 



Equation ([U describes our es timate of the mass l oss rate 
based on the formulation from Er kaev et al.l (120071) . Fxuv 
is the time- dependent total flu x between 1 - 1200 A, which 
IS given bv lRibas et all (l2005l) as a function of age for Sun- 
like stars, ^xuv is the planetary radius at which the at- 
mosphere becomes optically thick to X UV photons, which 
occur s at pressures around a nanobar dMurrav-Clav et all 
120091) . A^tide is a correction factor that accounts for the 
fact that mass only needs to reach the Hill radius to e scape 
(iLecavelier des Etangs etaH l2004t lErkaev et all 120071) . Fi- 
nally, e is an efficiency factor that parametrizes the fraction 
of the incident XUV flux that is converted into usable work. 
For this work we use e = 0.1 based on the observed photo- 
evaporation threshold described in Lopez et _a l. (2012). This 
value is similar to the efficiencies found byl Owen & Jackson! 
(l20T2h . Using more sophisticated photo-evaporation models 
for hot-Neptunes they found mass loss efficiencies that varied 
from 0.05 to 0.2. We make an additional conservative mod- 
eling choice by starting mass loss at an age of 10 Myr, since 
stellar XUV fluxes and planetary radii post-formation are not 
well understood at even earlier times. 

3. KEPLER-36: EXPLAINED BY MASS LOSS? 

Before we determine the possible mass loss histories of 
Kepler-36b & c, we must first examine their present day com- 
positions. For Kepler-36b, this is relatively straightforward 
since its density is consistent with a rocky compositions with- 
out a H/He envelope. In this case we find that it should be 
25.9±4 3 % iron, consistent with an Earth-like rock/iron ra- 
tio. This is derived by matching the observed mass and radius 
with our structure models in the absence of any H/He or wa- 
ter envelope. Likewise, the error bars include the observed Icr 
uncertainties in mass and radius. In contrast, Kepler-36c has 
much lower density and requires a substantial H/He envelope 
to explain its radius. To estimate its present day composition, 
we ran our thermal evolution models in the absence of any 
mass loss. To explain its current radius Kepler-36c must be 
8.6 ± 1.3% H/He, assuming an Earth-like core. To calculate 
the error bars we varied the planetary albedo from 0-0.8 and 
the heat capacity of the rocky core from 0.5-1 .0 J K"' g"' . We 
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Figure 1. Possible mass loss histories for Kepler-36b & c. The curves plot the mass fraction in the H/He envelope vs. time, while the colors indicate each planet's 
density at a given age. Likewise the inset shows the radius evolution of each planet. The black points on the right hand side indicate the current compositions of 
Kepler-36b & c, while the text lists the rock/iron mass predicted for each planet. Cun'ently, Kepler-36c requires ~ 8% of its mass in H/He, while Kepler-36b is 
consistent with an Earth-like composition. Nonetheless it is possible both planets fornied with ~ 22% H/He, but subsequent mass loss has distinguished them 
due to differences in their core masses. This provides a natural explanation for the large density contrast seen today in this system. Such evolutionary histories 
are a generic outcome of our evolution calculations, with a wide range of initial H/He masses and mass loss efficiencies. 



also included the observed uncertainties in mass, radius, cur- 
rent age, and incident flux. In addition, for Kepler-36b we can 
set an upper limit on the presence of any H/He envelope. It 
must be < 0. 1% H/ He, assuming a maximally iron rich core 
dMarcus et al.ll20Tol) . 

By knowing the present day compositions, we can deter- 
mine the mass of each planet's core. For Kepler-36b it is sim- 
ply the observed mass from TTVs, 4.5 ± 0.3 Mq, while for 
Kepler-36c it is 7.4±q 4 Mq. Using these core masses, we ran 
fully coupled models including both mass loss and thermal 
evolution, in order to determine the initial composition each 
planet had after formation. To ensure consistency, we check 
that these models with mass loss are still able to reproduce 
present day radius. 

Figure [T] shows the results of our coupled mass and thermal 
evolution models for both planets in the Kepler-36 system. 
The model assumes that both planets formed at their current 
orbits with the same initial H/He mass fraction. Each curve 
plots the fraction of each planet's mass in the H/He envelope 
vs. age, while the colors indicate the planet's density at that 
age. The black points at the right indicate the current age 
and composition of the Kepler-36 planets. We also show the 
radius evolution of both planets in an inset. Such large radii at 
young ages are a generic outcome of evolution models, since 
the H/He envelopes a re warm and are not degenerate. See 
iMordasini et al.l (120 1 2|) for a wider exploration of planetary 



radii for low-mass planets with H/He envelopes as a function 
of age. 

At 10 Myr, when we start photo-evaporation both planets 
are ~ 10 Rq. Kepler-36b rapidly contracts as it loses mass and 
is down to 3 by 100 Myr Kepler-36c is also vulnerable 
to mass loss. To have retained its current 8% of H/He Kepler- 
36c would need to have been formed with 22% H/He at 10 
Myr and a mass of 9.4 M®. Since Kepler-36b is consistent 
with no H/He envelope today, we can only set upper limits on 
its initial composition. However, if we assume that it formed 
with the same initial 22% H/He as 36c, then it would have lost 
its entire H/He envelope by the time it was 2 Gyr old. Thus 
the large discrepancy in the present day densities of the two 
planets, can naturally be explained by the fact that Kepler-36b 
is significantly more vulnerable to mass loss than Kepler-36c. 

However, this leads to the question, why is Kepler-36b so 
much more vulnerable than its neighbor? Although it does 
receive « 24% more incident radiation than 36c, this alone 
produces far too small an effe ct to explain the ne cessary dif- 
ference in mass loss histories ( Carter et al.ll2012l) . Instead, it 
is caused almost entirely by the difference in the masses of the 
two planets. As shown in L opez et al. ( 201Z) . planetary radius 
at fixed composition is relatively flat as a function of mass for 

1 -20 planets with significant H/He envelopes. As a 
result, density for these planets goes roughly like the mass. 
Since the mass loss rate is proportional to the inverse of the 
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Figure 2. This summarizes the resuhs of ~ 6000 model runs designed to explore how our coupled thermal and mass loss evolution models depend on basic 
planet properties. In each panel we have plotted models on a grid of incident bolometric flux, relative to the flux that the Earth receives from the Sun, vs. planetary 
core mass in Earth masses. At each point in the core mass-flux grid, we have over-plotted multiple models with different initial compositions. The size of each 
point indicates the initial H/He fraction a planet began with. Meanwhile, colors indicate how inuch of this initial H/He envelope is lost over time. Thus indigo 
and blue points are models where mass loss is unimportant while dark red points are models where the entire H/He envelope has been stripped off. Clearly, the 
mass loss history depends strongly on a planet's core mass as well as the incident flux it receives. In each case, there is a clear threshold region where mass loss is 
important but, at least some models are able retain significant H/He envelopes. This threshold region can be described roughly as a power-law with the threshold 
flux, Fth '^corc- Figure|4ji and equation l|6). The different panels show how the location of this threshold varies with age and mass loss efficiency. The top 
panels show the threshold at 50 Myr, 100 Myr, and 1 Gyr for our standard mass-loss efficiency of 10%. Most of the mass is lost in the first 100 Myr while alinost 
no mass is lost after 1 Gyr Likewise, the bottom panels show the results at 5 Gyr for eificiencies of 1.10, and 100%. Models with higher mass loss efficiencies 
are naturally more vulnerable to mass loss. 

possible that both planets started off with substantial amounts 
of H/He but that Kepler-36b lost its gaseous envelope to mass 
loss. 

4. THE ROLE OF CORE MASS 

4.1. A Parameter Study 

The Kepler-36 system clearly demonstrates that the mass 
loss evolution of planet depends on mor e than just the inci- 
dent XUV flux that a planet receives. In lLopez et al.l (120 12h 
we showed that there is a threshold in the observed popula- 
tion of transiting planets with measured densities and H/He 
envelopes. This threshold is well described by a critical mass 
loss timescale: 



of the average density 1/p, the mass loss timescale goes like 
MpX p orMp. 

However, rather than considering a planet's total mass to- 
day, it is useful to consider simply its heavy element mass 
today. Unlike the current total mass and bulk density, the 
heavy element mass should remain essentially constant as a 
planet contracts and loses its H/He envelope, providing a use- 
ful tool for estimating a planet's mass loss rate throughout 
its history. For simplicity, we assume here that the heavy el- 
ements are locked in a silicate/iron core, although it is also 
possible that some of these metals could be mixed into the 
H/He. Of course at a given composition, the mass of this 
rocky core Mcore is proportional to the total mass Mp. Thus 
the mass loss timescale f^i goes roughly like M^^^.^. We calcu- 
late that Kepler-36c's rocky core should be 65% more massive 
than that of Kepler-36b. As a result, we predict that the mass 
loss timescale for Kepler-36b was ^ 3 times shorter than for 
Kepler-36c. 

In Figure [1] we have presented one possible scenario for 
the evolution of Kepler-36. These results are insensitive to 
any reasonable variation in the mass loss efficiency or the ini- 
tial envelope fractions. For mass loss efficiencies significantly 
less than 0.1, it is not possible to construct a scenario where 
both planets b and c started out with identical envelope frac- 
tions. Nonetheless, even for lower efficiencies it is still quite 
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Here Fp is the incident bolometric flux that a planet receives 
from its parent star, is the current bolometric flux that the 
Earth receives from the Sun, and fxuv.Eioo = 504 erg s"' cm"^ 
is the XUV flux at the Earth when it was 100 Myr old. In 
iLopez et"an (120 12h we then performed a simple parameter 
study to show that this critical mass loss timescale was well 
reproduced by our coupled thermal and mass loss evolution 
models. 



In order to better understand in detail how mass loss evolu- 
tion depends on the mass of a planet's rocky core as well as the 
incident flux and mass l oss efficiency, we greatly expanded 
the parameter study from Lopez et al.l(l2012l) . Beginning at 10 
Myr, we ran over 6000 evolution models varying the incident 
flux, core mass, composition, and mass loss efficiency. We 
then recorded the radius, mass, and composition at different 
ages. The grid of initial conditions was evenly spaced loga- 
rithmically, with incident flux varying from 1 - 1000 F^, the 
rocky core mass from 1-64 M®, the initial H/He mass frac- 
tion from 0.1-80% H/He, and the mass loss efficiency from 
0.01-1.0. Figure |2] summarizes the results. In each panel we 
have plotted the incident flux vs. the core mass. Each circle 
corresponds to an individual planet that is color-coded by the 
fraction of its initial H/He envelope that is lost by the age in- 
dicated on the panel. At each point we overplotted multiple 
models with different initial compositions, in each case the 
size of the circle corresponds to the initial mass fraction in the 
H/He envelope. Finally the different panels compare results 
at different times and for different mass loss efficiencies. The 
top three panels show the results using our standard mass-loss 
efficiency at 50 Myr, 100 Myr, and 1 Gyr Meanwhile the bot- 
tom three panels show the results at 5 Gyr, a typical age for 
Kepler systems, for mass loss efficiencies of 0.01, 0.1, and 
1.0. 

Clearly, the mass loss history of a planet depends strongly 
on both the incident flux and the mass of the rock/iron core. 
The indigo models in the lower right of each panel have lost a 
negligible fraction of their initial H/He envelope. These mod- 
els have either relatively massive cores and/or receive little 
incident flux and so mass loss is unimportant to the evolution 
of planets in this part of parameter space. On the other hand, 
the dark red models in the upper left with low mass cores and 
high incident flux have completely lost their entire H/He en- 
velopes. This region of parameter space should be filled with 
highly irradiated rocky planets l ike CoR oT-7b and Kepler- 10b 
( lL6ger et al.ll2009l:IOueloz etal..2009-.Batalha et al.ll20TTl) . 

In between there is a transition region where mass loss is 
important but at least some of the models are able to retain 
a H/He envelope. This transition region is relatively narrow, 
spanning less than an order of magnitude in incident flux for 
a given core mass. For models in this transition region the 
relation between the size of the initial H/He envelope and the 
fraction of the envelope lost is extremely complicated and not 
always monotonic. There is trade off between the fact that 
planets with small initial envelopes have relatively little mass 
in those envelopes to lose, and the fact that planets with large 
initial envelopes have larger initial radii and therefore experi- 
ence much higher mass loss rates. For models at early times 
and/or low mass loss efficiency, the fraction of the envelope 
lost increases with the initial envelope fraction. On the other 
hand, at late times and/or high mass loss efficiencies it is the 
planets with small initial envelopes that are most vulnerable. 
As a result, much of the scatter in characterizing this tran- 
sition region is determined by variations in the initial H/He 
envelope fraction. 

In addition, the bottom three panels of Figure |2] make it 
clear that the location of this mass loss transition region de- 
pends on the mass loss efficiency. This makes sense; it is of 
course easier to remove more mass if the photo-evaporation 
is more efficient. Below, we will show that this behavior can 
also be well approximated as a power-law. In contrast, the 
behavior with age is more like an exponential decay. Most of 
the mass is lost in the first 100 Myr, while a comparison of the 
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Figure 3. The fraction of the initial H/He envelope which is lost to subse- 
quent photo-evaporation (/lost) vs. three different mass parameters, according 
to our models. Panel a) plots /iost against the mass of the rock/iron core while 
panel b) plots it against the mass of the initial H/He envelope and panel c) 
plots against the total initial mass. In each case, the different curves show the 
results for models with different initial H/He fractions /envelope > varying from 
1-60% H/He. All of these models receive 100 F^, assume e = 0. 1, start at 10 
Myr, and end at 5 Gyr Compared to envelope or total mass, core mass shows 
the least scatter between models with different composition. 

top right and bottom middle panels shows that relatively little 
mass is lost after 1 Gyr 

Likewise, in each panel the threshold for significant mass 
loss can be described as a power-law. If we define the cri- 
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terion for significant mass loss to be the evaporation of half 
the initial H/He envelope by several Gyr, then at a given age 
and mass loss efficiency this defines a power-law relation be- 
tween the rock/iron core mass and the threshold incident flux 
needed for substantial mass loss. The slope of this power-law 
is roughly Fth M^orc discussed in section |431 

4.2. Why Core Mass? 

The virtue of studying mass loss trends vs. core mass is that 
it allows us to define a narrow band of parameter space in Fig- 
ure|2]where the precise details of mass loss are important with 
relatively little scatter due to variations in the initial H/He en- 
velope mass. The currently observable planet properties, such 
as present day planet mass, radius, and density, are themselves 
highly dependent on the mass loss history and therefore on 
other unknown variables like the mass loss efficiency and ini- 
tial composition. This makes it difficult to separate the effects 
of variations in current planet mass from variations in other 
parameters like incident flux or mass loss efficiency. On the 
other hand core mass represents an initial condition which is 
constant throughout a planet's mass loss history. 

Rather than core mass we could alternatively chose to study 
trends against other theoretical parameters like the initial en- 
velope mass or initial total mass after formation. Like core 
mass, these parameters are independent of any subsequent 
mass loss evolution; however, in our models these parame- 
ters also less adept at predicting that mass loss evolution. In 
Figure [3] we show how /lost, the fraction of a planet's initial 
H/He envelope that is lost, depends on the core mass, the ini- 
tial envelope mass, and the total mass. Each of the different 
curves is for a different initial composition /envelope varying 
from 1-60% H/He. All the models are computed at 100 F^, 
assume e = 0.1, start at 10 Myr, and end at 5 Gyr 

Ideally we want a parameter that minimizes the scatter be- 
tween models with different initial compositions. Initial com- 
position is not directly observable, and determining a planet's 
mass loss history is a much more model-dependent exer- 
cise than determining its current composition. For observed 
planets like those in Kepler-36, it is possible to tightly con- 
strain the mass in a H/He envelope, even without the coupled 
thermal evolution models used here (iRogers & Sea ger 2010j). 
Moreover, as we show in Figures|4]|7] trends with composition 
are complicated and cannot be described by a simple power- 
law. Thus parameters that minimize the scatter between mod- 
els with different initial compositions are much better predic- 
tors of mass loss evolution. Comparing panels a) and b) in 
Figure |3] it is clear that plotting /ost vs. core mass produces 
vastly less scatter than plotting /ost vs. initial envelope mass. 
In panel c) we see that initial total mass is a reasonable pre- 
dictor of /lost, but with about twice as much scatter as in panel 
a). Unique among a planet's properties, M^oie is unchanged 
throughout a planet's evolution, is relatively model indepen- 
dent, and is a strong predictor of mass loss evolution. 

4.3. Scaling Relations for Coupled Mass Loss Evolution 

In order to quantify the dependency of mass loss on initial 
conditions, we examine how the location of the mass loss tran- 
sition region varies with each variable independently. This al- 
lows us to understand the qualitative behavior of our complete 
model in terms of a few simple scaling relations. Such rela- 
tions can be used for quick and relatively accurate estimates 
of the importance of mass loss for detected planets. 

We start by defining Fth as the threshold flux for a model 
with a given core mass, mass loss efficiency, age, and initial 
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Figure 4. Three panels showing how the threshold flux varies with core mass, 
mass loss efficiency and the initial composition. In panel a) we show how the 
threshold flux, the incident flux needed to remove half a planet's initial H/He 
envelope, varies with the core mass along with the best fit power-law. In panel 
b) we do the same for the mass loss efficiency. In panel c) we show how the 
threshold flux depends on the initial H/He envelope fraction. There is little 
dependence on composition until the planet starts out > 60% H/He. 
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Figure 5. Index of the core mass power-law from Figure |4]and equation (|6) 
vs. initial H/He envelope fraction. The dependence on core mass becomes 
steeper as the initial H/He fraction increases. 
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Figure 6. /lost the fraction of the initial H/He that is stripped after 5 Gyr 
vs. incident flux compared to Fth the threshold flux defined in equation j6). 
Below Fth. /lost increases roughly linearly with incident flux. Above ~ 2x 
fth. the envelop is completely stripped. 

composition to lose half of its initial H/He envelope. Holding 
all other variables fixed we examine how Fth varies with each 
parameter and try fitting a power-law relation. We then vary 
the other parameters across our entire parameter study and 
examine the scatter in these power-law fits. The dependence 
of Fth on Mcore, and /lost (the fraction of the initial H/He 
envelope that is lost) are all well fit by power-laws. As we 
previously described, the age dependence can be described as 
an exponential decay. For a 1 core and e = 0.1, the best fit 
exponential decay is: 

Fth exp (-(f - 140 Myr) /80 Myr)Fe + 3 AF^ (5) 

Meanwhile for systems older than 1 Gyr, age dependence is 
unimportant, allowing us to study Fth independently of age. 

Figure |4] shows in detail how the threshold flux and the 
amount of H/He lost depend on each variable in our parameter 
study according to our thermal and mass loss evolution model. 
For clarity we have picked representative values for the core 



Figure 7. f[„st according to the results of our full model from Figure |2]is 
divided by /lost according to the simple analytic desciiption in equatioiijSJ 
and plotted against initial H/He fraction. For initial compositions < 50% 
H/He, the two generally agree within a factor of 2 or better 

mass, mass loss efficiency, and the initial composition, but 
the results shown are generalizable across the entire parame- 
ter space. The default values correspond roughly to those for 
Kepler-36c: an 8 core, e = 0.1, and an initial composition 
of 20% H/He. Likewise, we choose models where 50% of a 
planet's initial envelope has been lost. We then vary Mcoie, 
and composition one at a time. For clarity, in each panel we 
list the variables that are being held constant. 

In Figure we examine how Fth depends on Mcore- This 
is well described by the over-plotted power-law with Fth oc 
^core- This power-law is closely related to the critical 
photo- evaporation timescale that we described in lLopez et alJ 
(l2012h . where the mass loss timescale f„,i goes like Mpp/Fp. 
This corresponds to the transition region in Figure |2] Because 
super-Earths and sub-Neptunes contain most of their mass in 
a rock/iron core, p, Mp, and Mcore all correlate strongly with 
each other, which implies that Fth oc Mpp oc M^ore- Our fitted 
power-law, Fth oc M^^^g, is slightly steeper than the simple 
dependence we would expect analytically. The difference is 
due to the slight dependence of radius on core mass. 

Likewise in Figure |4]3 we examine how Fth depends on e. 
This is also well fit by a power-law, however it is intrigu- 
ing that the dependence on e is not quite inversely linear as 
it would be in a simple mass loss timescale. This is one of the 
effects of coupling mass loss and thermal evolution. All other 
things being equal, as the mass loss efficiency increases, the 
mass loss timescale becomes shorter while the thermal cool- 
ing timescale is relatively constant. As a result, the radius 
decreases more slowly as mass is lost, meaning that slightly 
more mass will be lost over a planet's history. This means that 
/lost increases more than linearly with e. In turn this means 
that for fixed /ost, Fth decreases less than linearly with e. We 
can summarize both these trends in a single equation: 



flh = 0.5Fe 



Me, 



Mff 



2.4±0.4 



-0.7±0.1 



(6) 



Equation (|6]l describes the results of our power-law fits for 
the location of Fth. This is the incident flux a planet needs 
to receive from its parent star to remove half a planet's ini- 
tial H/He envelope over the planet's Hfetime, as a function of 
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Figure 8. Incident flux and core mass for 29 observed transiting planets with well defined masses less than 100 Mgj. Planets are color-coded by their current 
H/He envelope fraction. Red-brown open circles are consistent with being H/He free. The dashed line shows the Fth~W® relation from equations ^6] and ((T), 
scaled up slightly to account for complete stripping of H/He. Of the six planets that lie 1 cr to left of this relation, three are consistent with being rocky and three 
with having only water/steam envelopes. 



core mass Mcore, photo-evaporation efficiency e, and the flux 
that the Earth receives from the Sun F^. This was derived by 
fitting to all the models in our parameter study and the error 
bars correspond to the Icr errors in those fits. Unfortunately, 
the dependence on initial composition is significantly more 
complicated and cannot be described by a simple power-law. 

In Figure]?}; we plot Fth against the initial mass fraction in 
the H/He envelope at 10 Myr For modest initial H/He en- 
velopes, up to ~ 60% of the planet's total mass, increasing 
the H/He mass does not significantly affect fth- However, 
above this point the envelope's self gravity becomes suffi- 
ciently strong that density increases rapidly and it becomes 
increasingly difficult to remove any mass. At the same time, 
since these planets have such massive envelopes, removing 
a few Earth masses of H/He has a much smaller impact on 
their overall composition. Figure |5] shows the other main ef- 
fect of varying the initial composition. Here we show how the 
Mcore - Fth power-law index from Figure depends on the 
initial composition. In general, as we increase the initial en- 
velope mass, the radius-core mass relation becomes slightly 
steeper leading to a steepening of the Mcoi-g-Ff^ relation. 

In Figure |6] we complete the picture by showing how the 
amount of mass lost depends on all these aspects. We plot 
/lost, the fraction of the initial envelope that is stripped after 5 
Gyr, vs. the ratio of the incident flux and the threshold flux 
described by equation (|6]l. The relation is roughly linear as 
described by equation ([Tj, where if fi^^i > 1 then the planet is 



completely stripped. 



/lost = 0.5 



1.1 ±0.3 



(7) 



Finally, in Figure |7] we examine the quality of these simple 
analytic approximations compared to full results of our actual 
models. Here we plot the ratio of our full model and the an- 
alytic results of equations (|6]l and dTji vs. the initial composi- 
tion. The error bars represent the 1 a scatter due to variations 
in Mcore, e, and /ost- For initial compositions that are less 
then ^ 60% H/He, the two generally agree with each other to 
within a factor of 2. For more massive initial envelopes, the 
analytic description breaks down and overstates a planet's vul- 
nerability to mass loss. Also, it is important to keep in mind 
that these fitting equations are only rough approximations of 
the fully coupled evolution models and should not be used 
to make detailed predictions for individual planets. Nonethe- 
less, these equations are valuable in understanding the qualita- 
tive behavior of our complete model and in making statistical 
comparisons to large populations of planets. 

4.4. Comparison to Observed Population 

Figures 12] and |4] make a clear prediction about which plan- 
ets should be most vulnerable to mass loss. In addition to the 
incident XUV flux that a planet receives, its mass loss his- 
tory should depend strongly on the mass of its rock/iron core. 
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Planets that are either highly irradiated or have low mass cores 
are more vulnerable to losing any primordial H/He envelope. 
Thus we expect that we should not find planets with H/He 
envelopes above the Fth-A^core threshold relation in equation 
(|6]l. Any planets that are well above this relation should either 
be r ocky, or water wor lds which are less vulnerable to mass 
loss (iLopez et al.ll2012h . or have H/He envelopes so large that 
equations (|6]l and d?) breakdown. 

Figure |8] shows incident flux and core mass for all observed 
transiting planets with well defined masses less than 100 Mm 
and r adii less than 1.1 Rj from exoplanets.org (Wright etal] 
1201 ll) . Two planets were removed by the radius cut since they 
are known to be inflated and it is impossible to accurately de- 
termine a heavy element mass for these planets. To calculate 
core masses for these 29 planets we ran water-free thermal 
evolution models without mass loss for each planet. The er- 
ror bars on Mcore include the observed uncertainties on mass, 
radius, and age as well as theoretical uncertaintie s on the iron 
fracti on and thermal properties of the rocky core (ILopez et al.l 
120121) . Each planet is color-coded by its current H/He en- 
velope fraction. The red-brown open circles are planets that 
are consistent with being rocky today. The dashed black line 
shows the F±- Mcore relation from equation (|6]l scaled up by 
a factor of two so that equation (|7]i predicts complete strip- 
ping rather than only removing half the initial H/He envelope. 
Although the uncertainties are large, all the planets with sub- 
stantial H/He envelopes are consistent with being to the right 
of this threshold. Of the six planets that lie to the left of the 
threshold, three (Kepler- 10b, CoRoT-7b, and Kepler-20b) are 
consistent with being rocky and the other three with being 
water worlds (55 Cancri e, Kepler- 18b, Kepler-20c). 

5. EFFECTS ON PLANET RADII 
5.1. Trends in the Radius-Flux Distribution? 

Thus far we have examined the effects of thermal and mass 
loss evolution on planet mass and composition. However, for 
the vast majority on Kepler candidates neither of these can 
be determined. As a result it is worth examining the predic- 
tions that our parameter study makes for planet radii. Figure|9] 
shows the radius-flux distribution of the models in our param- 
eter study after 10 Gyr of thermal and mass loss evolution. As 
a reminder, we ran models on a log uniform grid of incident 
flux, core mass, and initial envelope fraction. Points show in- 
dividual models colored by their final H/He fractions while 
the grayscale boxes show the density of points in the distri- 
bution. In order to directly compare with other recent studies 
described below, we limited ourselves to core masses <15 
Mq, which also corresponds to th e high end of core masses 
that will trigger runaway accretion iMovshovitz et al.l (120101) . 
Modest variations of the maximum core mass, from 10-20 
M^, do not qualitatively change any of the results described 
below. Planets with more massive cores, will have likely un- 
dergone runaway accretion and have final radii > 4 R^,. In 
addition to varying the core mass, initial envelope fraction, 
and incident flux, we have also varied the iron fraction of the 
rocky core since this will smear out any trends with radius. 
We varied the iron fraction uniformly from pure silicate rock, 
to the maximum iron fraction allowe d by col lisional stripping, 
~60-80% in this mass range (Ma rcus et al.ll2010i) . For sim- 
plicity we restricted ourselves to only a single age and our 
standard mass loss efficiency 0. 1 . 

Our parameter study here is not meant to produce a realistic 
radius-flux distribution. Our choice of a log uniform distribu- 



tion in core mass, envelope fraction, and incident flux is meant 
to effectively probe the range of possible models. Nonethe- 
less there are key features in Figure |9]that should be observ- 
able in the flux-radius distribution of /Te/j/er candidates. First 
and foremost, there should be a decline in the rate of sub- 
Neptune sized planets, here defined as 1.8-4.0 Z?®, at high 
incident fluxes due to photo-evaporation. In particular, these 
planets should become comparatively rare for Fp > 100 Fq, 
which corresponds to perio ds <10 days. Such a n eff ect may 
have already been seen by lHoward etaH (120 12^) and'Youdin] 
(Hbl 1). At the same time there should be a corresponding in- 
crease in the frequency of rocky planets with < 1 .8 R^, due 
to the large number of stripped cores. 

lOwen & Wul ( 120131) recently performed a similar parame- 
ter study and predicted that there should be a significant gap 
in the radius distribution with few planets between 1.5 and 
2.5 /?0. The origin of this gap is easy to understand. Planets 
in this size range will typically be ^0.5% H/He. Such small 
envelopes are highly vulnerable to photo-evaporation; even 
tiny mass loss rates ^ 0.01 Gyr"' will be more than suf- 
ficient to strip these planets. Moreover, planets that formed 
with more substantial initial envelopes are unlikely to end up 
in this part of parameter space. This behavior was discussed 
at length in t he con text of evolution models for Kepler- 1 lb in 
ILopez et"an (120121) . If a planet experiences enough mass loss 
to remove several percent of its mass, then it is much more 
likely to lose its envelope completely than to end up with an 
envelope that is ~0.5% H/He. 

Nonet heless, we do not see a gap that is as clear cut as that 
found byl Owen '&Wiil(l20Tl . Instead we see a diagonal band 
in which models are relatively rare, although by no means ex- 
cluded. This "occurrence valley" is typically 0.5 R^ wide and 
occurs at slightly larger radii at higher incident fluxes. This is 
because in our parameter study we did not include planets that 
simply formed without any envelope at all. As a result, only 
planets that have lost their envelope to photo-evaporation end 
up being rocky. At the low flux end, only planets with initial 
envelopes ~0. 1 % H/He and the lowest core masses will lose 
their envelopes, resulting in relatively small stripped cores. 
On the other hand, at the high flux end even planets with ini- 
tial envelopes up to 1% H/He and core masses up to 10 
M0 are easily stripped. This removes somewhat larger H/He 
envelopes but results in larger stripped cores, moving the oc- 
currence valley up to ^2-2.5 R^. 

However even in the middle of the diagonal occurrence val- 
ley, we find models from our study. Partially, this is due to 
varying the iron fraction of the cores, which smears out radius 
trends by ^0.15 Rc^. Mostly however, it is due to running a 
comprehensive parame ter study that sample d a wide range of 
initial conditions. The lOwen & Wul ( 120131) study only tested 
five values for the core mass, without any initial composi- 
tions <1% H/He. In contrast. Figure |9]includes 200 different 
combinations of core mass and initial composition. In the ab- 
sence of photo-evaporation, this suite of models finely sam- 
ples the entire range of radii from 1-4 R^, as can be seen by 
the left most column in Figure |9] We suggest that the reason 
lOwen & Wu (2013) see a clear gap in the radius-flux distribu- 
tion is because their small sample of initial conditions do not 
adequately sample the parameter space. 

Thus far no such gap has been has been seen in th e ob- 
served distribution o f planet radii. iFressin et al.l (12013 ) and 
iPetigura et al.l (120131) recently performed careful studies of 
Kepler planet occurrence rates as a function of radius after 
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Figure 9. Final planet radius vs. incident flux in our parameter study. The greyscale boxes indicate the frequency of models that end in each box due to 
thermal and mass loss evolution; i.e., black boxes contain many models and white boxes few. At high incident flux there is a strong decrease in the frequency of 
~ 1.8 — 4.0 i?0 sub-Neptune sized planets and an increase in the frequency of < 1.8 rocky super-Earths. Points show the individual models, colored by their 
final H/He envelope fractions as result of photo-evaporation. The leftmost column of points at 1 closely approximates the distribution without any mass loss. 
Rust colored points in the bottom light coiTespond to rocky planets that have lost their envelopes. Just above these stripped cores there is a clear decrease in the 
frequency of models. Planets that enter this region have envelopes that are so small they tend to be stripped off completely. However, unlike iOwen & Wn<2013l) 
we do not see a clear gap in planet radius. 

Second, if sub-Neptunes typically form beyond the snow- 
line then in addition to rock, iron, and H/He, these planets 
could have large am ounts of water and other volatile ices 
jRogers et alj 1201 II) . Much like varying the iron fraction 
of the core, varying the water fraction could wash out any 
trends in radius, but to a much greater extent. A 5 planet 
that is 50% water will be ~0.5 .Rm larger than one with an 



correcting for false positives and the various selection effects. 
In both cases the studies find a flat occurrence rate below 
~2.8 with larger planets being significantly rarer Both 
these studies span a wide range of periods, out to 85 days for 
Ifressin et al. (2013) and 50 days for Petigura et al. (20131), 
and use fairly wide radius bins. As a result, it is perhaps 
not surprising that they would not detect our relatively nar- 
row occurrenc e valley. On the othe r ha nd, the flat occurrence 
rates found by ' Fressin et all (120131) and 'Petigura et al.' ('2013') 
seem i nconsistent with the wide gap proposed by Owen & Wu 

(Upo^. 

There are physical reasons why the occurrence valley seen 
in Figure |9]might not exist or be less pronounced. First, there 
could exist a large population of 1-10 rocky planets that 
simply formed without ever accreting a H/He envelope. This 
is possible if these planets formed through giant collisions af- 
ter the disk had already dissipated (Morbidelli et al. 20121). 
This would introduce another population of 1-1.8 plan- 
ets that would not show any strong dependence on incident 
flux due to photo-evaporation. Depending on how common 
these planets are, this could largely mask any gap in planet 
occurrence at radii <1.8 R^, which corresponds to incident 
flux <100 ^0 or periods longer than ~10 days. 



Earth-like composition dLopez et al.ll2012l) . Since our occur- 
rence valley is only ^0.5 R^ wide, varying the water con- 
tent sub-Neptunes from 0-50% would completely eliminate 
any dip in planet occurance. As a result, the presence or ab- 
sence of such a dip is a useful test for whether sub-Neptunes 
form in situ without an y water dChiang & LaughlinI 120131 
[Hansen & Murravll2012l) . or mig rate from beyond th e snow- 
line with large amounts of water ([Rogers et al.ll201 fl) . 

5.2. Relative Sizes in Multi-Planet Systems 

Beyond simply explaining individual systems like Kepler- 
36, models of mass loss evolution may shed light on many of 
the puzzles of planet occurrence statistics. One such puzzle 
is that many of the Kepler multi-planet systems like Kepler- 
11, Kepler- 18, & Kepler-36 exhibit regular ordering of their 
radii; i.e. each planet tends to be larger than the one interior to 
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it. This trend was recently quantified by ICiardi et"aLl (12013b 
which found that in pairs of planets from Kepler multi-planet 
systems there is a statistically significant tenden cy for the in- 
ner pl anet to be smaller than the outer planet. ICiardi et aP 
(120131) examined over 900 pairs of planets with periods rang- 
ing from 0.45 to 331 days and found that the inner planet is 
smaller in «60% of planet pairs. Furthermore, they found that 
the fraction of planet pairs where the inner planet is smaller 
rises to ^70% when the both planets are within 20 days, and 
that the trend is only apparent for planets that are smaller than 
Neptune. Assuming that these planets also have masses less 
than Neptune, these are precisely the planets that should be 
vulnerable to photo-evaporation. We find that this trend can 
be explained as a natural result of photo-evaporation. Unless 
there is a strong tendency for planets on shorter periods to 
have more massive rocky cores, then inner planets should be 
substantially more vulnerable to photo-evaporative m ass loss. 

On average we find that the inner planet in the Ci ardi et al.l 
(120 13h sample receives 7.9 times more incident flux than the 
outer planet. Moreover, since this trend exists for pairs of 
planets orbiting the same parent star, we know that they re- 
ceive the same XUV spectrum. Applying the scaling law de- 
rived in equation dTji, this implies that the inner planets in their 
sample should ty pically lose ^ 10 times as much H/He. In 
iLopez et"an ( 120 12l) we showed that there is a near one to one 
correspondence between radius and the H/He mass fraction. 
As a result, this increased vulnerability to mass loss should 
naturally lead interior planets to have smaller radii. In addi- 
tion, we would expect the fraction of pairs with smaller inner 
planets to rise at the shortest periods if these trends are in fact 
due to atmospheric mass loss. When both planet have periods 
shorter than ^^20 days, it is likely that both will be vulnerable 
to substantial mass loss and so there should be a larger impact 
on the relative radii. 

However, there are several factors that can diminish the im- 
pact of mass loss on relative radii. Firstly, some planets will 
be rocky super-Earths without any volatile envelopes. Like- 
wise our models predict that many planets on highly irradi- 
ated orbits should have envelopes completely stripped off. Of 
course, once a planet has lost its entire envelope its radius can 
not continue to shrink thus limiting any further differences in 
radii. Likewise, for planet pairs on less irradiated orbits, nei- 
ther planet might be vulnerable to photo-evaporation. Also as 
we have already shown, large differences in rocky core mass 
can overwhelm differences in incident flux. All of the effects 
combine to limit the usefulness of radii alone to understand 
differences due t o mass loss evolu tion and may explain why 
the trends seen by lCiardi et al.l (12013) are relatively weak. Un- 
fortunately, most of the planets in the Ciardi sample do not 
have mass measurements from radial velocity or TTVs, mak- 
ing it difficult to empirically test the importance of core mass 
on the trends they observe. 

Currently, there are 16 pairs of Kepler planets where both 
planets have well determined masses and meet the SNR and 
impact parameter thresholds described in Ciardi et al. (201 3). 
These include planets in Kepler-9 tHolman et al. 2010), 
Keple r- 10 (Batalha et al. 201 11), Kepler - 11 jLL ssauer et al. 
|2pir, 2013), Kepler-18 (Cochran et al. '2011^ Kepler-20 
( ^'essin et al. 201 2; Gautieretal. 2012), and Kepler-36 
(ICarter et al.ll2012l) . Of these 16 pairs, in four cases the in- 
ner planet is larger than the outer: Kepler-9b/c, Kepler- 1 1 c/f, 
Kepler- 1 Id/f, Kepler- 1 le/f. In all four of these cases the in- 
ner is significantly more massive than the outer planet. With 
only four cases the trend is not yet statistically significant. 



however, the tendency for inner planets to be either smaller or 
significantly more massive is a robust prediction of our mass 
loss models. Whenever the inner planet in a multi-planet sys- 
tem has a larger radius, it should also be significantly more 
massive. There should not be any highly irradiated pairs of 
planets where the inner planet is less massive than the outer 
but where the inner planet has enough of its mass in a H/He 
envelope such that its radius larger. 

It is worth mentioning that there are other processes in 
planet formation and evolution which could contribute to the 
ti-ends seen by Ciai'di et al. (2013). Ikoma & Hori (20121) 
showed that when super-Earths and sub-Neptunes form on 
highly irradiated orbits, the rate of H/He accretion is sig- 
nificantly slower when the local disk temperature is higher. 
Moreover, on short period orbits it is much easier for low 
mass planets to open a ga p in the disk, which wou ld also limit 
their envelope accretion (iHansen & Murrayl2012h . Assuming 
that planets in the Kepler multis formed in the same order- 
ing that they are in today, then inner planets should have had 
more difficulty in accreting large H/He envelopes. Likewise, 
given that proto-p lanetary disks evaporate from the inside out 
('Calv et et al.ll200 0). outer planets will have had more time to 
accrete an envelope. Nonetheless, given the critical mass loss 
timescale thr eshold that we iden tified observationally and the- 
oretically in 'Lopez et al.' (120 12i). and a concurring view ad- 
vanced by Wu & Lithwickl (120 12h . it seems quite reasonable 
that photo-evaporative mass loss plays an important role in 
producing the trends among planet pairs seen by Ciard i et aP 
(2013). Most likely, both planet formation and subsequent 
evolution combine to reduce the size of H/He envelopes for 
highly irradiated planets. 

6. SUMMARY 

There is growing evidence from both models and obser- 
vations that photo-evaporative mass loss plays an impor- 
tant role in the ev olution of highly ir r adiated super- Earths 
and sub-Neptunes dBaraffe et a l.l l2006t IJackson et al.l 120 12t 
IWu & LithwicS 120121) . iLecav elier Pes Etangs (2007) first 
proposed that there could be a critical mass loss timescale 
in the observed po pulation of hot Jupiters and Neptunes. In 
ILopez et"an (120121) we confirmed the existence this threshold 
for al l planets with meas ured densities down to 2 Mg. More- 
over, ILopez et"al] (120121) showed that this critical mass loss 
timescale is naturally reproduced by our coupled thermal and 
mass loss evolution models. Likewise, this mass loss thresh- 
old is also reproduced by other models which fully solve 
the hyd rodynamics of th e mass loss wind (Owen & Jacksoi^ 
2012; Owen & Wull2013h . H ere we have expand ed upon the 
parameter study performed in Lopez et"aLl (120121) and shown 
in detail how mass loss history depends on incident flux, core 
mass, and mass loss efficiency. 

We have shown that in addition to the amount of XUV 
irradiation received by a planet, the mass of its rock/iron 
core plays a critical role in determining a planet's photo- 
evaporation history. Moreover we have shown that this pro- 
vides a natural explanation for the large density contrast ob- 
served between Kepler-36 b&c. In order to better understand 
the role of core mass, we performed an extensive parameter 
study and provided approximate scaling relations which can 
be used for estimates of whether mass loss has been important 
for detected planets. 

Further, we showed that the compositions of the observed 
population of transiting planets are consistent with our de- 
tailed models and these scaling relations. Finally, we showed 
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that our coupled thermal and mass-loss evolution models 
make important predictions for the radius-flux distribution of 
Kepler candidates. In particular, we predict that sub-Neptune 
sized planets should become significantly less common at 
very short orbital periods. In addition, there may exist a nar- 
row "occurrence valley," which is a useful test for whether 
sub-Neptunes are formed in situ. We anticipate that future 
progress in this area will come from additional mass determi- 
nations of sub-Neptune size Kepler candidates, a better under- 
standing of XUV fluxes from all types of stars as a function of 
age, and further progress in modeling mass loss efficiencies in 
the framework of ID and 3D models. 
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